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The equations of motion describing buoyant fluids are often simplified using a set of approxima- 
tions proposed by J. Boussinesq one century ago. To resume, they consist in assuming constant 
fluid properties, incompressibility and conservation of calories during heat transport. Assuming 
fulfilment of the first requirement (constant fluid properties), we derive a set of 4 criteria for assess- 
ing the validity of the two other requirements in turbulent Rayleigh-Benard convection. The first 
criterion oA <C 1 simply results from the incompressibility condition in the thermal boundary layer 
(a and A are the thermal expansion coefficient and the temperature difference driving the flow). 
The 3 other criteria are proportional or quadratic with the density stratification or, equivalently 
with the temperature difference resulting from the adiabatic gradient across the cell A^. Numer- 
ical evaluations with air, water and cryogenic helium show that most laboratory experiments are 
free from such Boussinesq violation as long as the first criterion is fulfilled. In ultra high Rayleigh 
numbers (Ra > 10 16 ) experiments in He, one of the stratification criteria, scaling with a Ah, could 
be violated. This criterion garanties that pressure fluctuations have a negligible influence both on 
the density variation and on the heat transfer equation through compression/expansion cycles. Ex- 
trapolation to higher Ra suggests that strong violation of Boussinesq approximation could occur in 
atmospheric convection. 

PACS numbers: 47.27.Te Turbulent convective heat transfer , 44.25.+f Natural convection , 47.55.P- 
Buoyancy-driven flows; convection 

Keywords: Boussinesq approximation, Natural convection, Rayleigh-Benard, Turbulence 



I. INTRODUCTION 

The buoyant motion of fluids combines the non-linear 
complexity of hydrodynamics with additional degrees of 
freedom brought in by the thermodynamics. The exact 
equations for the buoyant flows are therefore very difficult 
to tackle directly. In 1903, Joseph Boussinesq introduced 
a set of approximations and derived simplified equations 
for buoyant motion [lj|. His three key approximations 
are : 

1) the fluid's properties are constant, 

2) the flow is incompressible (except for a linear temper- 
ature dependence of density yielding the buoyant term), 

3) the heat equation is not coupled to the flow mechanical 
energy ( "conservation of calories" ) 

Today, Boussinesq approximations are widely used to 
describe a large number of flows and to perform numer- 
ical simulations. A recent review presents an overview 
of works dedicated to the applicability of Boussinesq 
approximation and to its mathematical justification^. 
Oberbeck's name is occasionally associated with Boussi- 
nesq's to refer to the set of approximations presented 



above. As a side remark, it is worth noting that his 
convection model preserves compressibility and neglects 
all the time derivatives with respect to space derivatives 
Q.The resulting differential equations and physics sig- 
nificantly differs from the ones obtained by Boussinesq's. 

To be able to predict if Boussinesq's approximations 
are applicable, specific quantitative criteria are necessary 
for each type of flow. The present work derives such crite- 
ria for Rayleigh-Benard (RB) turbulent convection 0, [f|, 
assuming constant properties of the fluid (first approxi- 
mation satisfied). A vast literature is dedicated to this 
first approximation, in particular for high Prandtl num- 
bers in connection with mantle convection (for exam- 
ple see d, 0, HI, and references within). For the 
intermediate Prandtl numbers, recent progresses have 
been reported in the turbulent regime (for example see 
E3, EL IH OS and references within) . 

Boussinesq second approximation ("incompressibil- 
ity") and third one ("conservation of calories") consist 
in neglecting nearly half of the terms present in the ex- 
act equations of motion. When the magnitude of all the 
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terms can be estimated using typical scales of the flow, 
the applicability of the approximation can be predicted. 
In laminar RB convection, temperature and velocity gra- 
dients smoothly spread across the cell, which eases the 
choice of typical scales in the flow. A set of criteria for 
the applicability of Boussinesq approximation have been 
proposed in this low Rayleigh number (Ra) limit (for ex- 
ample see [l5l [l6j and ref. within). In turbulent RB 
convection (high Ra), choosing the correct typical scales 
is less straightforward due to the flow complexity. In 
the present work, we attempt to do so using various ex- 
perimental results and simple modelling. We derive a 
new set of 4 criteria for the validity of Boussinesq ap- 
proximation in a turbulent RB cell. In our derivation, 
we prioritise high Rayleigh numbers Ra and intermedi- 
ate Prandtl numbers (Pr) fluids (say Pr ~ 0.5 — 30), as 
this range of parameters is relevant to many geophysical 
flows of air and water, industrial flows and high Ra ex- 
periments with air, water, SF 6 and helium. Using these 
criteria, we investigate the validity of Boussinesq approx- 
imations in experimental set-ups with water, air and he- 
lium. In particular, we show that one criteria is much 
more stringent than the others in most experiments, at 
least below Ra ~ 10 16 typically. At higher Ra, a second 
criteria could become relevant. We also examine if de- 
viations from Boussinesq approximation can explain the 
apparent incompatibility of heat transfer measurements 
among very high Ra experiments. We find that it is not 
the case. 

The next section recalls the exact equations of mo- 
tion and singles out the different terms which have to 
be neglected to obtain Boussinesq's equations. In the 
third section, the magnitudes of these terms are com- 
pared with the retained terms by modelling turbulent 
convection. This yields a set of 4 criteria for the validity 
of Boussinesq approximation, which are summarised and 
discussed in the fourth section. The last section is a short 
conclusion. 



II. THE BOUSSINESQ EQUATIONS 

The temperature, pressure and density fields are de- 
composed into a constant (subscript "0"), a stratification 
(subscript "h") and a residual (9, p and Sp) contribution : 



T = T + T h + 9 
P = P +P h +p 
P = Pa + Ph + Sp 



(la) 
(lb) 
(1c) 



By definition, To is the average of the top and bot- 
tom plates' temperatures. Th is the linearised adiabatic 
gradient term defined by Eq. (f2"a|) . The pressure Pq + Ph 
and density po + ph fields satisfy the hydrostatic bal- 
ance Eq. (|2b[) and thermodynamic equilibrium Eq. (|2cp of 
a static stratified fluid with a temperature distribution 
T + T h (z). At cell's mid-height = 0), we set p h (0) = 



and P h (0) = 0. 

T h = 
dP h 



1 



dz 
dph 



ctgT 



-9 (Po + Ph) 



po + Ph dz 



t *Eb. 

dz 



dP h 

dz 



a 2 gT 



(2a) 
(2b) 

xg(po + ph) 

(2c) 



where a, g, c p and \ are t ne isobaric thermal expansion 
coefficient, the gravity acceleration, the isobaric heat ca- 
pacity per unit mass and the isothermal compressibility 
coefficient. 

The residual terms 9, p and dp account for fluctuations 
and for the inhomogeneities which are not captured by 
the stratification terms Th, Ph and ph- For example, 
fluctuations are generated by turbulence and inhomo- 
geneities result from the temperature forcing and confine- 
ment by the cell side walls. The potential temperature 
difference A and the adiabatic-gradient temperature dif- 
ference Ah across a convection cell of height h are defined 
as : 



A = 6{-h/2) - 9(h/2) 
A h = T h (h/2)-T h (-h/2) 



agToh 



(the notation V 2 will be used for the Laplacian operator 
to avoid confusion with the temperature difference A). 
The mass, momentum and energy conservation equations 



Vv 



1 Dp 

~p"dI 



p— = -VP + pg + F 



pc p 



Dt 
DT 

~Dt 



kV 2 T 



DP 

P* + *T— 



(4a) 
(4b) 
(4c) 



where D/Dt is the particle derivative and where F and 
$ account for the viscous force and corresponding dissi- 
pation and where k is the thermal conductivity. 

In this section, we assemble -step by step- a set of 
self consistent criteria which are sufficient to derive the 
Boussinesq equations from the conservation equations 
above. Since the successive approximations to derive 
Boussinesq equations are intricated, it is convenient to as- 
sume from the beginning that the density field is nearly 
homogeneous, or more precisely that the following two 
criteria are fulfilled (we omit the absolute values for bet- 
ter readability) : 



Ph < Po 
Sp < po 



(5a) 
(5b) 
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A. A preliminary thermodynamics relation 

Using the generalised Mayer relation and defining 7 = 
Cp/c v , we note that : 



1 



7 



PCpX Pogxh 



(6) 



This equation relates the density stratification result- 
ing from two different contributions : aA/, is the thermal 
expansion associated with the adiabatic-gradient and 
Pogxh is the compressibility associated with the hy- 
drostatic pressure. Since 7 is of order one or higher 
(critical point vicinity), this relation shows that the two 
contributions are comparable in strength. 



B. Momentum conservation 

Using the definitions Eq. (|lb[) and (fTc|) , the momentum 
conservation equation Eq. (|4b[) can be written : 



p—— = — V» ■ 
1 Dt F 



F - VP (l + (p + Ph)s 



By definition of Ph, the last two terms cancel out. An- 
ticipating on the solenoidal condition Vv = 0, the fric- 
tion is restricted to the shear term ??V 2 v which appears 
in Navier-Stokes equation (zero bulk viscosity). Using 
the criteria Eq.©, Boussinesq's momentum equation : 



Po 



Dt 



-Vp — poaBg + 77V v 



is recovered if the extra criterion is fulfilled 



Sp 



-p ab 



(7) 



(8) 



In the above equation, Sp is approximated at first or- 
der. We are assuming implicitly that higher order terms 
don't contain qualitative features which are significant 
for the flow dynamics. Such a "significant" feature is 
present -for example- in the viscous term (dissipation) 
and, consequently, it cannot be neglected with respect to 
the inertial term Dv/Dt, although the later is known to 
be larger in turbulent regimes. 



Mass conservation 



Using the definition Eq. lfTcj) . the mass conservation 
Eq. (|4a| can be written : 



Vv 



1 Dp _ 1 DSp _ v z dp h 
p Dt p Dt p dz 



Using the definition Eq. ([2c]) . the criteria Eq.© and 
IJ), and the thermodynamics relation Eq.©, we get : 



Vv 



DO v z p gx 



Dt 



Boussinesq's approximation requires that the velocity 
field is solenoidal (Vv = 0), which is satisfied if each 
term in the right-hand side (RHS) of the previous equa- 
tion can be neglected with respect to the typical veloc- 
ity derivatives dvi/dxi of the left-hand side (LHS) Vv, 
where i € {x, y, z} : 



7 dxi 
DO dvi 

Dt dxi 



D. Energy conservation 



(9a) 
(9b) 



Using the definitions Eq.© and ©, the energy con- 
servation equation Eq. p7i|) can be written : 

pCp TTt ~ kv2 ° = p ® + aT ^m + Vz9a [T ° P ~ T (po + Ph)] 

(10) 

Defining n = k/ pc p as the thermal diffusivity, Boussi- 
nesq's thermal equation (a "conservation of calories" 
equation) : 



D6 
~Dt 



- K v 2 6> = 



(11) 



is recovered if each term of the RHS of the equation 
Eq. (flO|) can be neglected. A sufficient condition for this 
is found using the criteria Eq. ([5]) and ([5]) and writing that 
To provides a reasonable order of magnitude of T. We 
find the new criteria : 



$/c p < 0' 
ocTq Dp 



poc p Dt 



v z ga 



[T h + 6 (1 + alb)] « 0' 



(12a) 
(12b) 

(12c) 



where e' e{ |f,,^, K 0} 



III. CRITERIA OF APPLICABILITY OF THE 
BOUSSINESQ'S APPROXIMATION 

In this section, we evaluate the criteria (O , © , ([5]) and 
(fT2| for the special case of high Ra (turbulent convection) 
and intermediate Pr. 

To check the first criterion Eq. (|5a|) , the equations (|2"c|) 
and (|2bp could easily be solved exactly but an order-of- 
magnitude solution for ph is enough for our purpose. Us- 
ing also Eq.©, we easily find that the criterion Eq. ([5a|) 
becomes : 



7 



< 1 



7-1 



< 1 



(13) 
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To estimate the other criteria, we discriminate between 
two regions of the cell where the flow phenomenologies 
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are known to differ : the bulk and the boundary layers 
near the top and bottom plates. Each criterion should be 
fulfilled in the two regions. Before proceeding with these 
two regions, it is useful to present order-of-magnitude fits 
of various quantities in the cell. 

A. Few orders of magnitude 

The scalings of various flow quantities have been re- 
ported in the literature of turbulent convection, in par- 
ticular in the Rayleigh-Benard cells. In this subsection, 
our aim is to present a few simple fits that will be used 
to estimate and compare the orders of magnitude of var- 
ious contributions. We are not suggesting here that the 
chosen power law fits have an underlying physical mean- 
ing. We neglect the inhomogeneities within the boundary 
layers and within the bulk of the flow, which result from 
finite size effects e.g. confinement of the large scale cir- 
culation. 

Below Ra = 10 12 and for 0.7 < Pr < 30, all reported 
Nu(Ra) data are in reasonable agreement (for studies of 
both Ra and Pr in this range, see for example [l7], [H, 
OS HI HI). We will use the fit : 

(Ra < 10 12 ) Nu ~ 0.065i?a 1/3 (14) 

The 1/3 exponent chosen in this fit is slightly larger than 
the exponents fitted in most studies (which are closer to 
0.30-0.32), even when the temperature drop in plates is 
compensated [22j |. But this simple fit still provides good 
numerical estimates for our purpose. 

Above Ra = 10 12 , the reported Nu(Ra) dependence 
fall into two sets : data from 2 experiments roughly fol- 
lows the scaling of Eq.fTTJ]) ([H,[2J|) while the others ex- 
hibit a steep increase of the Ra exponent from 0.31 up to 
typically 0.38 [II, [II, S3, The cnhanced-heat-transfer 
law can be fitted by : 

(10 12 <Ra< 10 14 ; Pr ~ 10) Nu ~ 0.015i?a 38 

The apparent incompatibility between the two sets of 
data is not understood and motivates present research 
(see also the numerical results [2jJ,|30|). Although all the 
data showing a heat transfer enhancement (including 
[3l| ) are consistent with the "ultimate regime" predicted 
by R. Kraichnan [32| . an extensive characterisation of 
this new regime is still missing. Therefore, we will not 
attempt to derive specific Boussinesq criteria for it in 
this paper. In the Discussion section, we will come back 
on this point to discuss very high Ra experiments. 

A large scale circulation or "wind" is known to occur in 
Rayleigh-Benard cells. The strength of the l.s.c velocity 
vi sc can be characterised by the fit : 

Re= hVlsc^ Q A5Ra l/2p r -3/4 (15) 
V 

which is a compromise between studies from various 
groups [H, [H, Hj, [H, [H, H3, HI] conducted with various 



operational definitions of Re for 0.7 < Pr < 30, based 
on local velocity or turn over time measurements [39j . 
An alternative estimate of the typical velocity in the 
bulk of the flow can be obtained from the rms velocity 
fluctuations. Comparing the fit of Re (above) with rms 
velocity data from the literature [H, H3, SO, EH we 
note that the scaling of the root-mean-squared (rms) 
velocity is approximately the same as the scaling of the 
large scale circulation velocity, and that the magnitudes 
are comparable (rms velocity is typically 40% of vi sc ). 
We will therefore consider only a single characteristic 
velocity at large scales. We should mention that, in 
the literature, the Ra dependence of Re has been fitted 
with effective exponents as low as 0.40 [35| (and even 
0.36 ± 0.05 at some location in a square cell [42|] ) an d 
as high as 0.65 (in the horizontal plane [37|). 

Many consistent studies report scalings of the rms tem- 
perature fluctuations in the bulk of the cell versus Ra 
(for example see @, H3, l43l lij. Hfl) but, to our knowl- 
edge, only Daya and Ecke [4fJ reports both a Ra and 
Pr dependence. The fit below is derived from Daya and 
Ecke original data. We should mention that the reported 
Pr dependence is poorly fitted by a power law but that 
such a fit still provides a fair order of magnitude for our 
purpose. 6* ms is defined as the rms temperature fluctua- 
tions normalised by the temperature difference A across 
the cell. 

0* ms ~ 0.3Ra-° 13 Pr-°- 5 (16) 

Finally, we will need an estimate of the time derivative 
d9/dt in the bulk of the flow. The references [H, 53] 
report corresponding measurements over a large range of 
Ra (5.10 s < Ra < 10 15 ) for Prandtl numbers of order 1- 
10 (see also [45[). They interpret a transition occurring 
near Ra ~ 10 11 with a change in power laws. Neverthe- 
less, within a factor 2, we can still re-fit all their data 
with the following power law, which is enough for our 
purpose (the Prandtl number dependence is unknown) : 

89/dt ~ Q.7Ra 0A ^- 

In the bulk of the flow, we expect advection to be the 
most effective heat transport mechanism, which implies 
88/ dt ~ vV#. The relation shows that we can also es- 
timate the time derivative from i>V£?, taking v ~ vRejh 
and using recent temperature-gradient measurements for 
10 9 < Ra < 10 10 and Pr ~ 5.4 |, Averaging the 
measured gradient from the various locations inside the 
cell|48|. we obtain the fit : 

80/ dt ~ vV0 ~ t;[V0| =s 70 Ra- a 1165 Re^- 

nr 

It is satisfactory that both estimations are in reasonable 
agreement (magnitude and Ra exponent). We will favour 
the first fit because its has been validated over a wider 
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range of Ra. It is convenient to define a dimensionless 
derivative 9* dt as : 



'dt 



89 h 2 n , 

i9i Ak 



(17) 



where a; = (Pr/5) x accounts for the unknown Prandtl 
dependence through the unknown parameter X. For in- 
termediate Pr, we will assume that x is of order 1. 



B. Boundary layer region 

The present knowledge of scaling in the boundary 
layers is not as extended as in the bulk, in particular 
due to the complex interplay between the thermal and 
the velocity boundary layers (for recent progresses, see 
H3, [EH)- Nevertheless, it is well established that 
each boundary layer concentrate nearly half of the to- 
tal temperature drop A over a typical distance of order 
A = h/2Nu. An internal velocity scale vu can be ob- 
tained writing that total heat flux experiences a cross- 
over within the thermal boundary layer from a molec- 
ular diffusion transport (right on the plate) to a con- 
vective transport at a distance A from the plate, which 
gives Vbi = ft/A = 2Nu.n/h. An external velocity scale 
is also imposed on the thermal boundary layers by the 
large scale circulation present in the bulk of the cell 
vi sc = Re.v/h. Using the estimates of Nu and Re given 
previously, we note that vi sc > vu- In the following, we 
will choose the velocity scale which gives the more strin- 
gent criteria ("sufficient conditions"). Using these tem- 
perature, length and velocity scales, we evaluate below 
the criterion ([5b). ©, © and ■ 

Local thermodynamic equilibrium implies the differen- 
tial equation : 

dp = —padT + pxdP 

Subtracting the stratification contribution (f2"c)) from 
this local equilibrium, the criteria Eq. (l5b|) and © are 
found to be equivalent to : 



a9 <C 1 
XP -C ct9 



(18a) 
(18b) 



In the boundary layer, we have 9 ~ A/2. The order 
of magnitude of p is assessed taking the divergence of 
Eq. ([7]) and using the solenoidal property of the velocity 
field. We find a Poisson equation for the pressure : 

09 

V 2 p=-poV{{v\7)v}+ap g— (19a) 

8z 

The first term of the RHS gives the usual dynamical pres- 
sure of order : 

PO f 2 , 2n PO 2 

p ~ y( u &* + y isc ) ~ y VIsc 



and the second RHS term, a buoyant source of pressure 
variation, can be estimated by integration as : 



V ■ 



gp aAX 



Substituting the expressions vi sc = vRejh and A = 
h/2Nu and using the fits for Nu and Re, we find that the 
buoyant contribution to the pressure is smaller than the 
kinematic contribution. Neglecting the former and using 
Eq.©, both criteria Eq. (fT5)) becomes, up to a factor 2 : 



«A«1 
PrRe 2 aA h 



Ra 1 - 7- 1 



< 1 



(20a) 
(20b) 



The second criterion Eq. (|20b[) is comparable or less 
stringent than Eq. (|13p . so we can drop it. 



In the boundary layer, the next criterion Eq.jy.ii is 
estimated as : 



7 A 



After simplification by the typical velocity v and using 
the expression of A, we find that the above condition 
is less stringent than the criterion Eq. (fT3|) . We will 
therefore ignore it too. 



To evaluate the criterion Eq. (|9b[) in the boundary 
layer, the particle derivative of the temperature is es- 
timated with a diffusion process (heat is indeed trans- 
ferred between the plates and the fluid by diffusion in 
the boundary layers) : 



D9 
~Dt 



A/2 
AV^ 



(21) 



and the space derivative of the velocity is estimated 
from the typical internal velocity as Vbi/X. The external 
velocity vi sc is not chosen because it would give a more 
loose condition. Substituting the expressions of D9/Dt, 
vu and A in the criterion Eq. (|9b[) gives criterion Eq. (|20a[) 
again. 

We now consider the three criteria Eq. (Tl2")) associated 
with the energy equation. The three RHS term 9' of 
these inequalities can be estimated with the character- 
istic scales used above. All turn out to have the same 
typical magnitude, as expected from the chosen bound- 
ary layer model : 



e' 



89 
dt 



89 

'8x4 



8^1 

dxf 



kA 

2X 2 



We estimate the LHS of the inequalities Eq. (JT3J) using 
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vi sc as the typical velocity ( "sufficient condition"). 



< 



— ( Vlsc \ ' 
c„ V A / 



aTp Dp aT p vi sc 2 /2 
p Cp Dt p Cp \/vi sc 



v z ga 



T h + 6(1 + aT )] 



< 



vi sc ga 



Ah 
2 



A 



(1 + aT ) 



Substituting the expressions of vi sc (Re) and X(Nu) in 
the previous equations, we find up to prefactors of order 
unity : 



A h 2Pr 2 Re 2 

x 5 <C 1 

Jo na 

Pr 2 Re 3 
aA h x — - < 1 



Ah (Ah 
T { A 



=?■ =^ + vT Q + 1 



2RaNu 
PrRe 



4Nu 2 



< 1 



(22a) 
(22b) 
(22c) 



The second condition (|22b[) is comparable or more 
stringent than Eq. (fT3")) and will replace it. The aTp term 
in the parenthesis of Eq. (|22c|) gives a condition com- 
parable or less severe than (|22b[) and can therefore be 
dropped. The unity term in the parenthesis gives a con- 
dition that can be compared to Eq. (|22ap and it turns out 
to be typically up to 10 times more stringent for low Ra 
(Ra ~ 10 7 ), comparable for Ra ~ 10 13 , and less strin- 
gent above. We will come back on these 2 conditions 
later. The Ah/ A term of Eq. (|22cp gives the condition : 



A? PrRe 
T n A ANu 2 



C. Bulk region 



(23) 



In the bulk, for the same reasons as in the boundary 
layers, the criteria Eq. (|5bf and JH]) are equivalent to : 



a$ < 1 

xv *c 

The first condition is weaker than the criteria Eq. (|20ap 
since the magnitude of temperature fluctuations in the 
bulk is significantly smaller than the total temperature 
drop across the cell (9* ms <1). 

To assess the second condition, the magnitude of 
the pressure p is estimated using the Poisson equation 
Eq. tfTTj)) . The buoyant contribution is estimated by anal- 
ogy with electrostatic, apogdO / dz being the pressure 
source with a negligible remote contribution from the 
boundary layers [63] and a local contribution from the 
bulk. 



< 



Po 



apog\ 



,86 



where I is a correlation length of the temperature gradi- 
ent 89 /dz in the bulk. The buoyant contribution |§||^ 2 
can be estimated from the temperature rms fluctuation 
A9* ms , the typical temperature derivative 88/dzi, and 
using the inequality 

t d8 9 ,89 ,89 AO* 

Estimating the typical velocity as v = vRe/h and using 
the fit of Re and 9* ms , we find that the buoyant contri- 
bution to the pressure variation is comparable or smaller 
than the kinematic term for Ra > 10 s . We will ignore 
it. A sufficient condition for the criteria \P a 9 in the 
bulk is therefore : 



1 



7" 



PrRe 2 
2RaO* 



< 1 



(24) 



This new criterion is comparable to Eq. (|22b|) . We 
choose the new one because it is derived from tighter 
estimates. 

To estimate the criteria Eq. d^aj) . we take v z ~ vRe/h 
as the typical vertical velocity. The typical derivative 
dvi/dxi is estimated assuming that the velocity fluctua- 
tions in the bulk of the flow are driven by the large scale 
circulation on a large scale h and follows a Kolmogorov- 
like cascade. Then, the energy injected into cascade is re- 
lated to the viscous dissipation $, which is proportional 
to the square of velocity spatial derivatives ( 



>roportion; 
BE El) : 



JfL = Re* x — ~ $ ~ 
h ft 4 

from which we get the estimate 



8x t 



Re 3 ' 2 — 
h 2 



Using those estimates and Eq.©, it is straightforward 
to show that Eq. ([5aj) in the bulk is less stringent than 
Eq. tfTB"]) . We can therefore ignore this condition. As 
a comment, the LHS of Eq. (|5a|) has been estimated 
using the magnitude of the velocity at large scales 
while the RHS has been estimated at small scales where 
the derivative is known to be larger. The solenoidal 
condition Eq. ([5aj) is nevertheless still accounted for at 
large scale only, as can be checked estimating the RHS 
as vi sc /h (one gets Eq. (fT3|) ). 



To estimate the next criterion Eo. (|9bj) . we need an 
upper bound for the particle derivative D9/Dt. We can 
consider the relaxation by molecular diffusion of a A-thick 
plume of excess temperature A/2. The result is the same 
as in the boundary layer : Eq. (|2"Tj) . The criterion Eq. (|9bD 
becomes in the bulk : 



2Nu 2 
PrRe 3 ' 2 



aA« 1 



Using typical fits Re(Ra), Nu(Ra), we find that the 
previous inequality is less severe than (|20a[) and can be 
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ignored. 

We now consider the three criteria Eq. (fT2| associated 
with the energy equation. Boussinesq thermal equation 
Eq. (fTTj) states that the advection of temperature inho- 
mogeneties is only affected by the molecular thermal dif- 
fusion. In the boundary layers, we showed that these two 
transport processes have the same efficiency. In the bulk, 
we expect advection to become the dominant process, ex- 
cept in some localised regions (on the thermals interface) 
where diffusion could be as efficient as in the boundary 
layers. To estimate the three inequalities Eq. (fT2]) , we 
will require that each LHS term is significantly smaller 
than the advection term vV# ~ d6/dt = O^An/h 2 of 
the RHS. It should be noted that the molecular diffu- 
sion term in the thermal equation provides a cut-off (or 
"dissipation" ) mechanism for temperature fluctuations at 
small scales and that we are not aiming to compare the 
magnitude of this cut-off term with the LHS terms of 
Eq.®. 

To estimate the LHS of the first thermal criterion 
Eq. (|I 2a[) . Eq. (|2"5)) provides an order of magnitude of 
the viscous dissipation $ which is found consistent with 
the Direct Numerical Simulation (DNS) results of Verz- 
icco [4lj. The criterion Eq.fT2aJ in the bulk becomes 
Eq.||2BaJ). 

To estimate the next criterion Eq. (|12b[ ). we assume 
that pressure and temperature fluctuates with similar 
typical time scales. The derivatives Dp/Dt and 86/ dt 
of Eo. (|12b[) are replaced with their integrated value p ~ 
p v 2 /2 and 6 ~ A9* ms . We find Eq.JUQ. 

To estimate the LHS of the third thermal criterion 
Eq-KUSp, we take v z ~ Re.v/h, 9 ~ A6* ms and T h ~ 
Ah/2 to find Eq.(|2"oc]). 



A h Pr 2 Re 3 



To 
aA h x 



RaO* m 

PrRe 2 
2Ra9* a 



Ah 
To 



A h 



2A6± 



aT + 1 



PrReOt 



'dt 



< 1 



•C 1 



< 1 



(26a) 
(26b) 
(26c) 



Using the typical orders of magnitude of Re, 
and 9g f , we now compare these three thermal conditions 



AO* 

rms 



Eq. (|26|) with previous ones. The first condition is less 
stringent than Eq. (122cp (unity term) at low Ra and more 
stringent than Eq. (|22ap at high Ra. We choose to drop 
the condition Eq. (|22a|) and to merge Eq. (|22c[) (unity term 
only) and Eq. (|26a[) into : 



A h /PrRe 
— x - 

T \ANu 2 



Pr 2 Re 2 
RaO 



at 



< 1 



(27) 



Up to a numerical factor (1 — 7 _1 ) of order one, the 
second thermal condition Eq. (|26b[) is equivalent to the 
criterion (|24|) . We now consider one by one each term in 
the brackets of Eq. (|26cp . The unity term gives a condi- 
tion which is less stringent than Eq. (|2"T)) . The aTo term 



gives a condition less stringent than Eq. (|26b[) and can be 
ignored. The A h /A9* ms term gives a condition compa- 
rable, or more stringent for high Ra, than Eq . (|23| . We 
therefore keep the former one and re-write it alone as 



PrRe 



To A 29 



< 1 



(28) 



dt 



D. Summary of the criteria and discussion 

The set of the most stringent criteria for the applica- 
bility of Boussinesq approximation leaves us with only 
four sufficient conditions Eq.([20ajl. ([24]) . (|27]) and ([28]) . 
which are recalled below, up to a factor 2. 



Ah 
To 



and we recall that 



aA h x 
Pr 2 Re 3 

To A 



aA < 1 
PrRe 2 

<C I 



Ra9* ms 

PrRe\ 

ANu 2 J 

PrRe 
x 

2e dt 



< 1 



< 1 



(29a) 
(29b) 

(29c) 

(29d) 



aA h = a 2 gT h/c p = poffxM 1 ~ 7 l ) ~ Po9X h 

These 4 criteria are a first outcome of this paper. 
Although the choices of these specific criteria is not free 
from some arbitrariness, the order of magnitude of the 
LHS of these inequalities should be fairly robust to these 
choices. 

Physically the first criterion is required for incom- 
pressibility in the boundary layer. The other criteria 
are more difficult to relate to a specific physical effect. 
For example, the complete heat equation has a term 
accounting for the thermo-mechanical ener gy coupling 
which directly results from stratification [64j |. This 
term -which appears on the LHS of Eq. (|12cjl - gives 
contributions proportional to aA/,, Ah/Tq and A$/TqA. 
The fourth criterion above accounts for the A^/TqA 
contribution alone. The third criterion above includes 
the Ah /To contribution but also another condition 
necessary to neglect the viscous dissipation with respect 
to advected heat flux. The second criterion above allows 
to neglect variations of density due to the pressure 
variations, versus those due to temperature variations 
fEq. (|2~4")) . It also allows to neglect the cooling/heating 
associated with the pressure variations experienced by a 
fluid particle in the heat transport equation (Eq. (l22bp 
and Eq.j26b|). 

These 4 criteria don't allow to assess a-priori the fulfil- 
ment of the Boussinesq approximation for given Ra and 
Pr since it requires some knowledge of Re, 9* dt and 9* ms . 
A further step in the analysis consists in replacing these 
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quantities by fits versus Ra and Pr. For example, us- 
ing the fits presented previously, the 4 previous criteria 
becomes : 

aA « 1 (30a) 



O.l-Ra '" (aA h ) < 1 



Ra 01 Pr-°- 25 x 
200 



10i?a- Ui 'Pr 



0.17 D 0.25 



0.1Ea 01 Pr°- 25 z ( ^ ] « 1 



r 

A? 



(30b) 
<Cl (30c) 



(30d) 

where we rounded prefactors and exponents to recall 
that the above inequalities are derived from order-of- 
magnitude fits. We recall that the order of magnitude 
of x is taken as unity for Pr of order 1-10. 

These equations highlight parameters having a direct 
influence on the validity of the Boussinesq approxima- 
tion. For example, the dependence versus Ra is very 
weak. The equations Eq. ([30|) indicates that high Ra 
are weakly unfavourable but the uncertainty in the fits' 
exponents prevents from drawing precise conclusions on 
this weak dependence. For given Ra and Pr, we found 
general features. Tall cells are unfavourable because of 
stratification (linear and quadratic dependence versus h 
coming from Ah). Close proximity to the critical point 
is unfavourable too since a has a strong divergence in a 
region where Ah stays approximately constant. Finally, 
large temperature differences A have opposite influence 
on the criterion Eq. (|30ap (unfavourable) and on the 
criterion Eq. (|30d[) (favourable). In practice all these 
parameters can't all be varied independently. 

Numerical applications are now proposed for a few 
typical laboratory Rayleigh-Benard cells. Let's first con- 
sider a water cell of height h = 1 m, mean temperature 
To = 50° C and driven with A = 30° C (near atmospheric 
pressure). Such cells are being studied -for example- in 
physics labs in Lyon, Hong-Kong and Santa Barbara 
(for example see [54|, [55|, [56[). We found that the criteria 
Eq. ([3T)]) arc all fulfilled. The most stringent one is by far 
Eq. (|30ap , which is nevertheless still valid with 2 decades 
separation between the two sides of the inequality. The 
three other criteria are fulfilled with at least 5 decades 
of separation. 

In h = 8 m air cell at atmospheric pressure, such as 
the one operated in Ilmenau (for example see (57j). we 
also find that the first criterion Eq. (|30a[) is the most 
stringent. For A = 50°C and T = 50°C, we find 
a ration 6 separation between the two sides of this 
inequality. The other criteria are fulfilled with 2.5 or 
more decades of separation. 

Although atmospheric convection differs from RB 
convection of dry air, it is interesting the examine 



the criteria Ea.(|30|l for h = 10 km, A = 50 C and 
T = 0°C Using dry air properties, we find that 
Eq. (|30b[) is violated by more than one decade and than 
the LHS of Eq. (|30cp is close to 1, indicating that the 
Boussinesq approximation is not valid in this academic 
case. Surely these numbers should be taken with much 
caution because they rely on the extrapolation up to 
Ra ~ 2.10 21 of fits derived at much lower Ra and 
probably in a different turbulent regime. 

The case of Helium cells, such as those being operated 
in Grenoble and Trieste cannot be addressed with the 
same general procedure as above because, from one run 
to another, the temperature and the fluid density are var- 
ied and the parameter a does change accordingly (while 
Ah/h remains close to 2.10 -3 - 3. 10 -3 Km^ 1 ). Besides, 
the modelling presented in this paper ignores the possi- 
bility of occurrence of the Kraichnan's ultimate regime 
of convection [32j , while the interpretation of some high 
Ra data is done according to this regime. These words of 
caution written, we can still get guidelines from the de- 
rived criteria. We numerically find that the two first cri- 
teria Eq. (|30a[) and Eq. (|30b[) are the most stringent ones, 
while the two others are fulfilled by at least 3 orders of 
magnitude for cells of height 20 cm and 1 m Therefore, 
we focus below on the validity of the second criterion for 
experimental data fulfilling the first one. 

Data from \24 . [2 5l 26 1 have been corrected from the 
side wall effect [58l l59l. l6Cll| and filtered to reject those 
with aA > 0.2, that is a 10% maximum relative variation 
of the density in each boundary layer (first criterion). All 
cells are cylindrical and their height and aspect ratio Y = 
diameter / height are reported in the legend of figure[T] 

The top plot shows the LHS of the criterion Eq. (|30b[) 
versus Ra. We recall here that the fit of 8* ms used to 
derive Eq. (|30b[) from Eq. (|29b|) is smooth over nearly 6 
decades, up to nearly Ra = 10 15 and that we assume that 
its extrapolation over one extra decade (up to Ra = 10 16 ) 
provides a correct order of magnitude. From the plot, we 
find that this non-Boussinesq parameter reaches values 
larger than 0.2 for the highest Ra (Ra ~ 10 16 ), which 
may no longuer be considered as significantly smaller 
than unity. This suggests that Boussinesq approxima- 
tion could be violated for such ultra high Ra. 

Besides, this shows that increasing the He cell height h 
to reach higher Ra -in Boussinesq conditions- is at best 
efficient like h 2 and not h 3 . Indeed, if we neglect the 
weak Ra dependence of Eq. (|30b|l and increase by a fac- 
tor -say- 10 the cell height, a should be 10 times smaller 
for aAh to stay constant (fulfilment of Eq. (|30bp \ On the 
other hand, A can be made 10 times larger (fulfilment of 
Eq. (|29a[) ). In temperature-density area where these high 
Ra data are obtained, a 10 fold decrease of a is roughly 
correlated with a 10 fold increase of the thermal diffusiv- 
ity but with little change in the kinematic viscosity. All 
parameters considered, the maximum "Boussinesq" Ra 
will only be 100 times larger. One the other hand, this 
implies that such tall cells could be interesting models 
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E. Conclusion and Perspectives 





10° r 




10- 1 r 


< 




a 




CO 


10" 2 : 


.1 Ra° 


l0- 3 r 


o 


I 

10" 4 r 



C Chavanne el al. [1997] (1=0.5, h=20cm) 
• Niemela et al. [2000] (1 =0.5. h=100cm) 
. T Niemela el al. ]2003] (1=1, h=50cm) 








1 








Ra 







10"' 



CO 






ra 






DC 


10^ 




Z3 






Z 








We derived a set of 4 criteria (Eq.(|2"9"|0 for assessing 
the applicability of Boussincsq approximation in a fluid 
with constant properties at high Rayleigh number Ra 
and intermediate Pr. The fourth criterion, scaling 
quadratically with stratification, was found to be always 
less stringent than the 3 others in our numerical estima- 
tions. 

The first criterion qA <C 1 simply reflects the in- 
comprcssibility condition in the boundary layers. The 
3 others are indirectly associated with stratification 
effects. They are written as functions of the adiabatic 
gradient but the thermodynamics relation Eq.([6]) allows 
to translate them into equivalent density stratification 
criteria. Contrary to a widespread idea, tall cells are not 
necessarily more Boussinesq than small ones for given 
Ra and Pr, because of these stratification effects. 

Using fits derived from the RB literature, we found 
that the first criterion aA <§; 1 is the most relevant one 
for assessing the "Boussinesqness" of typical laboratory 
experiments conducted with water and air. In geophys- 
ical flows for which stratification occurs on very large 
scales, the second and third criteria could become more 
relevant than the first one. 



FIG. 1: Data from three cryogenic He experiments [24j, |2g, l2g , 
restricted to those fulfilling aA < 0.2. The top plot shows the 
parameter appearing on the left of the Boussinesq criterion 
Eq. ()30b[ l. The Boussinesq condition is better fulfilled for the 
values of the parameter smaller than one (the Chavanne et al. 
data near Ra ~ 10 12 sitting below 10 -4 have been obtained 
in liquid helium). The bottom plot shows the compensated 
Nusselt NuRa~ 1/3 versus Ra. We note that, for Ra ~ I0 13 
(vertical line), the datasets already sits on different Nu(Ra) 
trends while some of their non Boussinesq parameters, dis- 
played on the top plot, roughly overlap. 



for atmospheric convection. Indeed, according to the es- 
timate presented above we expect a similar violation of 
the Boussinesq approximation in the atmosphere. 

The bottom plot of this figure illustrates that the 
compensated Nusselt NuRa^ 1 / 3, has different trends for 
Ra > 10 12 as mentioned earlier in the text. Regard- 
ing this issue, the top plot also shows that the non- 
Boussinesq parameter is below 10 -2 for Ra ~ 10 13 and 
-more important- that it is approximately the same for 
datasets with different Nu(Ra) trends (the vertical line 
across the plots illustrates this point). This shows that 
non-Boussinesq deviations of the sort encompassed by 
the present analysis cannot explain the discrepancy of 
heat transfer efficiency among He cells for I0 12 < Ra < 

io 14 m, h m, he m m, nil. 



In cryogenic helium, in addition to the first criterion, 
the validity of second one (Eq. (|29b[) ) deserves to be care- 
fully examined at ultra high Ra, above 10 16 typically. 
This point and its consequences on very high Ra exper- 
iments is discussed in the text. At such high Ra, this 
second criterion is a sufficient condition for density fluc- 
tuations in the bulk to be more correlated to temperature 
fluctuations than to pressure fluctuations (Eq. lfM)) ) and 
to neglect, in the heat equation, a thermal heat source 
(or sink) term due to pressure-fluctuation-induced com- 
pression (or expansion) (Eq. (|26b[ )). 

At lower Rayleigh numbers (typically IO 12 < Ra < 
I0 14 ), we find that the apparent incompatibility of heat 
transfer measurements reported in literature cannot be 
explained by such non-Boussinesq effects. 

In the Boussinesq approximation, the equations of the 
fluid are only parametrised by Ra and Pr. The present 
work shows that a weaker approximation of the equations 
-preserving some terms proportional to the parameters 
a A and/or a Ah (or equivalently XP9^)~ would enable to 
account for first deviations from the Boussinesq limit of 
turbulent Rayleigh-Benard cells filled with a fluid with 
constant properties. As a perspective, a DNS simulation 
conducted with such generalised equations at moderate 
turbulent Rayleigh numbers (say IO 8 < Ra < IO 10 ) but 
with large enough values of the non-Boussinesq parame- 
ters aA or aAft would provide valuable quantitative in- 
formation on deviations from the Boussinesq limit. In 
particular, simulation with aAj, of order unity can still 
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be performed assuming incompressibility. Such simula- 
tions would give a precise quantitative meaning on the 
<C symbol used in the proposed Boussinesq criteria. 
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